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ABSTRACT 

We re-examine the role of collisions in decoupling the H I 21-cni spin temperature 
from the cosmic microwave background (CMB). The cross section for de-exciting the 
21-cm trasition in collisions with free electrons is more than 10 times larger than it is 
in collisions with other atoms. If the fraction of free electrons in the diffuse cosmic gas 
is between 10 and 30 per cent then collisions alone can decouple the spin temperature 
from the cosmic microwave background (CMB), even in moderately under-dense re- 
gions at z > 15. This decoupling is especially important during the very early stages 
of re-ionization when a Lya continuum background had yet to be established. As a 
detailed example, we develop a semi-analytic model to quantify 21-cm emission signa- 
tures from a diffuse gas which is partially ionized at z > 10 by an X-ray background. 
We find 21-cm differential brightness temperature fluctuations with a mean of > 8 mK 
and a rms value as large as 5 mK, for a frequency resolution bandwidth of 100 KHz 
and a beamsize of 3 arcmin. Another example where free electron-atom collisions are 
important is during the recombination of bubbles ionized by short-lived UV sources. 
When the ionized fraction in these bubbles drops to 10-20 per cent their differential 
temperature can be as high as 10 mK. 

Key words: cosmology: theory - intergalactic medium -large-scale structure of the 
universe 



1 INTRODUCTION 

Recent observational data of the universe have substantially tightened our grip on the cosmological model and its fundamental 
parameters. Most prominent of these data are maps of the cosmic microwave background (CMB), in particular those measured 
by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite (e.g. Spergel et al. 2003), the 2dF and SDSS galaxy redshift 
surveys (Percival et al. 2002; Zehavi et al. 2002), and the Lya forest seen in spectra of quasars (QSOs) (Croft et al. 1999, 
2002; McDonald et al. 2000 & 2004, Nusser & Haehnelt 1999, 2000; Kim et al. 2004; Viel, Weller & Haehnelt 2004). Despite 
the remarkable recent achievements, our understanding of the initial stages of the formation of structure remains poor. The 
weakest link is the era of hydrogen re-ionization at high redshifts. The WMAP polarization measurement (e.g. Spergel et al. 
2003) imply a universe that must have been ionized at redshifts z > 13. However, the imprinted signal of re-ionization in the 
CMB is insensitive to the details of the initial stages. A probe of these details can be found in H I absorption features in spectra 
of high redshift QSO. These features can reveal the spatial distribution of ionized regions as a function of redshift (Becker et 
al. 2001, Fan et al. 2002). Unfortunately, bright QSOs are hard to detect at high redshifts and, also, the inference of constraints 
on the three-dimensional structure of re-ionization can be very tricky (e.g. Nusser et al. 2002). A more promising probe of 
reionization, especially of the initial stages, is observations of the redshifted 21-cm emission/absorption lines produced by H I 
in the high redshift universe (Field 1959, Sunyaev & Zel'dovich 1975, Hogan & Rees 1979, Subramanian & Padmanabhan 
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1993). The 21-cm line is produced in the transition between the triplet and singlet sublevels of the hyperfinestructure of 
the ground level of neutral hydrogen atoms. This wavelength corresponds to a frequency of 1420MHZ and a temperature of 
T« = 0.068 K. If ni and no are the populations of the triplet and singlet ground state of HI atoms then the spin temperature 
is defined by ni/no = 3exp(— T./Ts). The CMB radiation tends to bring Ts to its own temperature of T^-^^ Ki 2.73(1 + a) K 
(Mather et al. 1994). An H I region would be visible against the CMB either in absorption if Ts < T^gj^ or emission if 
Ts > T^BR • Two mechanisms for decoupling the spin and CMB temperature have been proposed (Field 1958). The first 
is the Wouthuysen-Field process (Wouthuysen 1952, Field 1958) which is also termed "Lya pumping". In this process the 
hyperfine sublevels are interchanged by the absorption and emission of a Lya photon. The final product of this process is 
a spin temperature that equals the "color" temperature of the Lya photons. The color temperature eventually approaches 
the thermal temperature of the gas as a result of recoil effects by repeated scattering off the H I atoms (Field 1959). The 
other mechanism is spin exchange by means collisions of atoms with free electrons and other atoms. Both of these effects 
work towards bringing the spin temperature close to the gas thermal temperature*. The efficiency of Lya pumping has been 
studied thoroughly by several authors (Scott & Rees 1990, Madau, Meiksen & Rees 1997, Ciardi & Madau 2003, Tozzi et al. 
2000, Chen & Miraldar-Escude 2004). For reionization by stars alone the Lya pumping is believed to be very efficient even 
when the volume filling factor of ionized regions is loss than a few per cents. The reasoning behind that is as follows. During 
the early stages of reioinization all UV photons emitted with hv > 13.6 ev by a stellar source are consumed in ionizing the 
surrounding gas (mainly H I). A Lya continuum photon with 10.2 < hv < 13.6 ev is free to travel until it is redshifted to a 
hf « 10.2 ev where it is captured in the Lya line of a ground state H I atom. Such a photon travels a comoving distance of 
Ra ~ 6000 h-^Mpc(l + z)"^''^(l-(4(/w'/13.6)/3)^/^) (about 350 h^^Mpc at z = 20 for hv = 13.6 ev) before its capture. Since 
a star emits at least 4 times more photons in the Lya continuum than in the UV (e.g. Ciardi & Madau 2003)), one expects 
a significant amount of the continuum photons to reach the H I far beyond ionized regions. According to Ciardi & Madau 
(2003) Lya pumping should already have become efficient when the volume filling factor of ionized regions is as small as a few 
per cent. This conclusion is invalid when the ionizing sources axe highly clustered. The Lya continuum from an aggregation of 
sources could decouple the spin-CMB temperature in a sphere of maximum radius of Ra around it, independent of the number 
of Lya continuum photons emitted by these sources. However, the ionized region around this aggregation is proportional to 
the number of sources in it. Further, X-ray photons from accreting black holes could easily have played an important role in 
an early re- ionization epoch (e.g. Ricotti & Ostriker 2004). Because of the large uncertainties in describing reionization (c.f. 
Cen 2003, Wyithe & Loeb 2003) it is prudent to examine in detail the consequences of other mechanisms for decoupling the 
spin-CMB temperatures. 

In the absence of free electrons, spin exchange by H I atom-atom collisions is important only in regions with high density 
contrasts, S > 20[(1 + z)/10]"^ (e.g. Ciardi & Madau 2003). This condition is satisfied in minihalos with Tcsi, < Tvir < lO"* K 
which could contain a significant fraction of H I at « ~ 8 (Iliev et al. 2003) in the absence of external heating sources. In this 
paper we argue that collisions of H I atoms with free electrons and other atoms in moderate density regions with 70-90 per cent 
neutral fraction can indxicc significant 21-cm emission against the CMB. As an example of this wc focus on pre-reionization 
(i.e. early stages of reionization) by X-ray photons (Oh 2001, Venkatesan, Giroux & ShuU 2001, Machacck, Bryan & Abel 2003, 
Ostriker & Gnedin 1996, Ricotti & Ostriker 2004, Ricotti, Ostriker & Gnedin 2004). X-ray photons have such a small cross 
section for H I ionization so that they rapidly form a uniform ionizing background (e.g. Ricotti & Ostriker 2004) . The average 
ionized fraction in the presence of this background is in the range 5-20 per cent for 30 > « > 10, which is consistent with 
the current constraints on the observed soft X-ray background (Dijkstra et al. 2004). Under these circumstances, collisions 
are very efficient at decoupling the spin-CMB tempaxetures. We develop a simple semi-analytic model to estimates the 21-cm 
emission fiuctuations induced by spin-CMB decoupling by means of collisions in the framework of the X-ray scenario. We 
also consider 21-cm emission during the recombination of "bubbles" ionized by sporadic short-lived sources of UV ionizing 
radiation. Gas in a newly formed ionized bubble would have a temperature of a few 10'' K and would immediately start cooling 
(mainly by Compton cooling) and recombining into H I. Electron-atom collisions inside the bubble maintain the recombined 
H I at Ts S> Tf^BR- This high Ts amounts to a significant 21-cm emission against the CMB when the H I fraction is the range 
70 — 90 per cent. This emission could last for ~ 10*Yr so that repetitive formation of these bubbles could yield an appreciable 
signal. 

The outline of the paper is as follows. In §2 we summarize the basic equations. In §3 we present results for the 21-cm 
differential brightness temperature at mean gas density as a function of temperature, the H I fraction, and redshift. The 
temperature fluctuations axe estimated numerically in §4. With the exception of subsection 4.3, all of §4 is devoted to the role 
of collisions in the X-ray pre-reionization scenario. We discuss the results and conclude in §5. 



* There is actually another decoupling mechanism. Immediately after hydrogen atoms recombine they populate the triplet and singlet 
hyperfine sublevels according to the ratio 3 : 1, leading to Ts = oo. However, this decoupUng is just a transient as very rapidly the CMB 
radiation field brings Ts to finite values. 
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2 THE BASICS 

Wo assume throughout that spin exchange by means of colhsions is the dominant mechanism for decoupUng the spin temper- 
ature from the CMB. In this case The spin temperature is given by (Field 1959), 

rj. _ ^cBR + .'A-^K , . 

1 + yc ' 

where 

Cio T* 

yc = ^7f^, (2) 

^10 J K 

with Cio denoting the total de-excitation rate as a result of collisions with free electrons and other atoms. Following Field 
(1959), we write 

2/c =yHl'^m +J'e"e • (3) 

In the following we adopt the values of y^j and listed in Table 2 of Field (1958). The corresponding values for Cio are in 
good agreement with other published data for the relevant temperature range considered here (Allison & Dalgarno 1969). At 
Tk = 3000 K and lO"* K the ratios of to are 20 and 13.8, respectively. 

Intensities, /(i^) at radio frequency are expressed in terms of brightness temperature defined as Tb = I(v)c^ /2kv^ , where 
c is the speed of light and k is Boltzman's constant (Wild 1952). The differential brightness temperature (hereafter, DBT) 
against the CMB of a small patch of gas with Tg at redshift z is (e.g. Ciardi & Madau 2003), 

l + z\ /0.3~ 



5Ti 



16mKa;jjj(l + 5) [l 



TcBn\ f^^^h\ [f^ + ^\ ^0.3^] ^ 



/ V0.02 



w J \n 



where S and the gas density contrast and the fraction of H I in the patch, respectively. 



3 THE BRIGHTNESS TEMPERATURE AT MEAN DENSITY 

We first study the DBT, JTb, at mean density [5 = 0) as a function of redshift. We assume throughout that the spin and CMB 
temperatures are decoupled only by means of electron-atom and atom-atom collisions. In Fig. 1 we show a contour map of 5Tb 
at mean gas density (5 = 0) at a = 15 as a function of x^^ and the gas temperature Tk- The DBT peaks at x^^ « 0.72 and 
Tk ^ 3000 K. It is insensitive to Tk above 3000K, but varies appreciably as a function of x^j. For 10" K < T K < 3 X 10* K 
there is an increase by a nearly factor of 2 between Xjjj = 1 and 0.7. 

In the top panel of Fig. 2 we examine the spin temperature, Ts , as a function of the neutral fraction, x-^^ . This temperature 
is very sensitive to x^j because of the enhanced coupling between Ts and Tk produced by electron-atom collisions. It is well 
below Tk, but significantly larger than the CMB temperature T^b^ = 2.73(1 -|- z). The middle and bottom panels in Fig. 2 
show (5Tb as a function of for Tk ~ 10* K and 10^ K, respectively. Here also we see that a small free electron fraction can 
substantially boost up (5Tb. For x^^ « 1 collisions with free electrons become less effective at decoupling Tg from T^^^, making 
5Tb oc (1 — Tpgj^/Ts) (see Eq. 4) drop by a factor of 2 fi-om its peak value at Xjjj = 0.8. Note that the (5Tb curves peak at 
ajjjj — 0.7 — 0.9, close to the values expected in the X-ray pre-reionization scenario (Ricotti & Ostriker 2004). In Fig. 3 we show 
the redshift evolution of (5Tb at (5 = for x^^ = 0.8 and 0.9, as indicated in the figure. In a realistic scenario for X-ray pre- 
reionization varies with redshift (Ricotti & Ostriker 2004). The history of X-ray pre-ionization, if indeed happened, could 
in principle be probed by observations of the mean temperature as a function of observed frequency, f — 1420MHz/ (1 + z). 
In practice, this signal maybe overwhelmed by foreground source (Shaver et al. 1999; Baltz, Gnedin &i Silk 1998, Ricotti, 
Ostriker &c Gnedin 2004). 



4 THE 21-CM TEMPERATURE FLUCTUATIONS 
4.1 The model 

As a concrete example for assessing the role of collisions we consider the DBT in the X-ray pre-reionization scenario (Oh 2001, 
Venkatesan, Giroux & ShuU 2001, Machacek, Bryan & Abel 2003, Ostriker & Gnedin 1996, Ricotti & Ostriker 2004, Ricotti, 
Ostriker & Gnedin 2004) . The mean free path of X-ray photons is much larger than the correlation length of the sources and 
their mean separation. Therefore, in this scenario X-ray photons act as a uniform ionizing background and the H I fraction, 
Xjjj, at a point with density contrast, 5, is determined by the equation for photoionization equilibrium, 

rx^,=a^i\n{l + S){l-x^,f , (5) 
© 0000 RAS, MNRAS 000, 000-000 
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Figure 1. Contour map of the DBT as a function of the neutral fraction, a^pjn ^'^'i ^^"^ thermal temperature, Tk- The calculation 
is for a gas at mean density at redshift z = 15. The DBT levels (in mK) are indicated by the contour labels. 



where V is the X-ray photoionization rate per hydrogen atom, \ is the recombination cross section to the second excited 
atomic level, and n is the mean density of total (ionized plus neutral) hydrogen. According to Ricotti & Ostriker (2004), the 
mean H I fraction should be maintained at a level Xjjj ~ 0.8 — 0.9 starting from redshifts 20 — 40 in order to match the WMAP 
value of Te ~ 0.17 for optical depth to Thomson scattering (Spergel et al. 2003). Therefore, here we fix the ration of^^n/T by 
assuming that is given at mean density, i.e. at 5 = 0. Then x^^ for any 5 can be easily found using (5). 

The DBT fluctuations over the sky and as a function of frequency contain significant cosmological information. In the 
case of the X-ray scenario, these fluctuations are a direct probe of the density fluctuations during the pre-reionization period 
since, as can be seen from Fig. 4, the DBT is an increasing function of the local density. In the following we make an 
attempt at assessing the statistical properties of these fluctuations. Ideally one would like to have a suite of hydro simulations 
combined with radiative transfer in order to compute the ionized fraction as a function of redshift and position for a given 
cosmological model (e.g. Gnedin & Shaver 2004, Ricotti, Ostriker & Gnedin 2004). Unfortunately, current simulations of 
this type are limited to small boxes (~ 10 — 20 h~^Mpc) making them unsuitable for detailed statistical analyses. Here we 
adopt the alternative approach of modeling the dark and baryonic matter using semi-analytic approximations. We work in 
the framework of the ACDM cosmological model with fim = 0.3 and S^a = 0.7. In order to generate a prediction for the 21-cm 
emission maps in this model we follow the following steps, 



(i) generate a random gaussian realization in a cubic grid box. This field is to represent the dark matter (DM) density field 
in a ACDM cosmology. 

(ii) obtain the gas density field by filtering the DM density with a smoothing window designed to mimic pressure effects in 
the linear regime. 

(iii) assume that the gravitationally evolved (non-linear) gas density field is related to the linear field by a log-normal 
mapping so that the probability distribution function of the evolved field is log-normal (Bi, Boerner & Chu 1991). 

(iv) compute the ionized fraction at each point in the box for a uniform X-ray background. The amplitude of the ionizing 
background is fixed by the value of the ionized fraction at mean density. 

(v) compute the spin temperature at each point assuming a uniform gas temperature of 10'* K (e.g. Ricotti & Ostriker 
2004) throughout the box. 

We work with a cubic grid of 256"^ points in a box of size L = 40 h~^Mpc (comoving) on the side. In k-space the smoothing 
window we use to mimic pressure effects is of the form (Bi, Boerner & Chu 1991), 

1 + (fc/fej) 
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Figure 3. The DBT at mean density as a function of redshift, plotted for two values of Xjjj as indicated in the figure. The gas temperature 
is 10* K. 

In this expression kj is the Jeans scale which depends on the gas temperature and redshift as, 



where Cg is the speed of sound. We smooth the DM density using the filter (6) with fcj corresponding to Tk = 10* K in order 
to derive the linear gas density contrast, 5l- The log-normal model then gives the non-linear gas density, 5, as (Coles & Jones 

1991, Bi 1993), 



where (Tl is the rms value of Jl- This transformation yields a log-normal probability distribution function for 5. The value of 
ctl determines the amplitude of the evolved gas density. Here we tune ctl so that the desired rms value, as, oi 5 is attained. 
There is a large uncertainty in the determination of os and, also, in the form of the filtering window used to infer the gas 
density from the DM distribution (e.g. Bi & Davidsen 1997, Hui & Gnedin 1997, Nusser 2000, Viel et al. 2002). In any case, 
one of our goals is to show that a detection of the 21-cm DBT fluctuations could actually constrain the properties of the gas 
density field as long as 05 is high enough. We will present results for as = 1 and as — 0.6. These values are inferred from 
the non-linear power spectrum of the ACDM computed from the linear power spectrum using the recipe outlined in Peacock 
(1999). The linear power spectrum is normalized such that the linear rms value of of density fluctuation in spheres of radius 
8 h~^Mpc is ag — 1. The non-linear power spectrum is smoothed with the filter given in (6) to derive an estimate for as- 
The outcome of this is as close to 0.7. However, because of the uncertainties in the shape of the filter and the value of as 
this value of ag can only be a a crude estimate. Hence we present results for as = 1 and 0.6. Once the evolved density field is 
obtained on the cubic grid, we solve the equations (1), (4), and (5) to derive the DBT, JTb, at each grid point. 

4.2 The fluctuations 

We begin with analyzing two dimensional maps of the DBT. The beamsize, A^, and the frequency bandwidth, Au, deter- 
mine the resolution at which "an observed map" represents the underlying temperature fluctuations. We ignore foreground 
contamination and refer the reader to Oh & Mack (2003) and Di Matteo et al. (2002) for details. We also do not include 
the effect of redshift distortions in the analysis (e.g. Nusser 2004). Here we mimic observations of DBT as follows. We place 
the box at the desired redshift and arbitrarily pick the plane xi — X2 in the box to represent the plane of the sky in a given 
direction. Using the angular diameter distance, we express the coordinates xi and X2 as angular positions in the sky. The X3 
coordinate is expressed in terms of the frequency 1420MHz/[l -|- z{x3)], where z{x3) is the redshift corresponding to 0:3. The 
DBT in the box is then given as a function of the angular position (in the xi — xi plane) and frequency (the xz direction). 




(7) 



1 5 = exp (5l - cr£/2) , 



(8) 
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Figure 4. The DBT as a function of tlic gas density contrast at z = 15. The H I fraction at 5 = is Xjjj = 0.8. For 5^0 we compute 
ajjjj assuming a uniform ionizing background as explained in the text. 

At redshift 15 our box {L = 40 h~^Mpc) represents a region spanning a 2.59MHz range in frequency and subtending 19.22 
arcmin on the sky. We smooth the DBT field with an eUipsoidal gaussian window having of a FWHM of Av in the frequency 
axis, and a FWHM of in the perpendicular (angular) plane. 

Figs 5 and 6 arc contour maps of the DBT in the plane of the sky for gas rms density fluctuations of cr^ = 1 and as — 0.6, 
respectively. All maps in all panels of the two figures show the DBT in the same slice of the box for a beamsize of Ad = 3.4 
arcmin. In each figure, the panels in the top row are for Af = 10 KHz, while the bottom row is for A,^ = 100 KHz. The panels 
in the column to the left are for neutral fraction (at mean density) of x^^ — 0.8, while the column to the right shows maps for 
ajjjj = 0.9. As a quantitative guide to the map we list, in table (1), the mean, < 6Th >, and the rms, < 6T^ >^^^. A visual 
comparison between the top {Au = 10 KHz) and bottom {Av = 100 KHz) panels in the figures (5) and (6) reveal substantial 
differences. The fluctuations in the top panels arc significantly higher and appear to show more structure. The comoving 
physical scale of a A6 = 3.4 arcmin is about 6.5 h^'^Mpc, nearly 3 times larger than the physical scale corresponding to 
Av = 100 KHz. Despite the large beamsize, small scale fluctuations are clearly seen for Av = 10 KHz. 

There are clear differences also between the contour maps for = 0.9 (panels to the right) and for = 0.8. These 
arc caused by the electron-atom collisions which are more effective for x^^ = 0.8. The higher fluctuations amplitude seen in 

the panels to the left is consistent with figures (1) and (2). 

The amplitude of the gas density fiuctuations clearly plays an important role in determining the level of DBT fluctuations. 
According to table (1) the DBT fluctuations for cr^ = 1 are about 50% higher than in the ag = 0.6 case (see also figures (5) 
and (6)). However, the mean values of STh are smaller in the former case. The reason is the larger volume of space occupied 

by under-dense regions for as = 1. 

In figure (7) we plot the rms of the DBT, < ST^ >^^^, as a function of the beamsize, A9. The curves correspond to various 
cases as explained in the figure caption. The rms values as a function of are closely related to the average of angular 
correlation function within an angular separation of ~ A6 (e.g. Peebles 1980). None of the curves reveal any characteristic 
feature related to a particular physical scale. This is not surprising since no such scale can be identified in the X-ray pre- 
reionization scenario employed here. With the exception of the cases corresponding to Au = 100 KHz with as = 0.6 (the 
solid and dashed lines with the square symbols), all curves show a non- negligible level of fiuctuations. However, an inspection 
of the contour maps in figures (5) and (6) shows that the DBT fluctuations between some regions can significantly be larger 
that the rms values. 

We turn now to the fluctuations along the frequency axis. These are of interest since the correlation along the frequency 
axis can be extremely useful for disentangling the cosmological signal from foreground contamination (Gnedin & Shaver 2004, 
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Table 1. The mean and rms value of the DBT for Xjjj = 0.8 and 0.9, as computed from the gas density fields with ag = 1 and as = 0.6. 
In each entry the first number in the bracket is the mean < (5Tb > and the second is the rms < >^/^ . 



^Hi ~ 0-8 = 0.9 

0-4 = 1 (9.05, 8.94) (7.86, 7.44) 
as = 0.6 (9.29,5.60) (8.05, 4.59) 




Figure 5. Contour map of the fluctuations in the DBT estimated from the gas density field with cr^ = 1 at z=15. The panels in the left 
and right columnn, respectively, correspond to Xjjj = 0.8 and 0.9 at (5 = 0. The top and bottom panels are for frequency bandwidth of 
A.v = 10 KHz and 100 KHz, respectively. The maps are for an angular resolution of 3.4 arcmin and the box size is 40 h~^Mpc on the 
side. The solid thick line marks the 8K level of the DBT. Fluctuations below and above the 8K level are represented by the light solid 
and dotted lines, respectively. The contour spacing is 5 mK. 
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Figure 6. The same as the previous figure but for the gas density field with = 0.6. The contour spacing in this map is 3 mK. 



Zaldarriaga, Furlanetto & Hernquist 2004). In the top panel of Fig. 8 we plot the DBT fluctuations in a segment of 40 h~^Mpc 
through the numerical box. All curves show DBT in the same segment for a bandwidth of 100 KHz, but for different angular 
resolutions. As is clear from the bottom panel, the DBT is significantly correlated up to separations corresponding to a few 
fi-^Mpc. 



4.3 Brightness temperature in recombination bubbles 

In this paper we have focused on 21-cm emission induced by collisions during a pre-reionization period in which X-rays are the 
main source of ionizing photons. An alternative scenario would be pre-reionization by UV radiation produced by short-lived 
sources (of lifetime much shorter than the recombination time-scale). 

In this scenario re-ionization begins with the production of individual ionized bubbles engulfing the sources. While a bubble 
recombines electron-atom collisions decouple the spin and CMB temperature efficiently so that significant 21cm emission is 
possible. The sources also produce Lya photons which can decouple the CMB from the spin temperature of H I. But, if the 
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Figure 7. The rms value of the DBT. The circles are attached to curves obtained from the density field with erg = 1, while the squares 
are for the fields with ag = 0.6. Solid, dotted, dashed and dash-dotted lines, respectively, correspond to (Ai/, a;jjj) = (100 KHz, 0.8), 
(10,0.8), (100,0.9), and (10,0.9). 



filling factor of the bubbles is small and the birth rate of these short lived sources is low then Lya photons may not become 
important at the early stages of re-ionization. In a separate paper we will present a detailed analysis of a situation like where 
the source are identified with massive first stars (e.g. Omukai & Nishi 1998, Abel, Bryan & Norman 2000, Bromm, Coppi & 
Larson 2002) . Here we only consider the differential antenna temperature in a region undergoing recombinations after it has 
been completely ionized by a short-lived source. 

To obtain brightness temperature as a bubble recombines we have numerically integrated the recombination equations 
assuming Compton cooling (e.g. Peebles 1993) to be the dominant cooler at at the relevant redshifts. The recombination 
coefficients as a function of temperature are taken from Storey & Hummer (1995). The de-excitation coefficients for atom- 
atom and atom-electron collisions are from Field (1958), which are in good agreement with other calculations in the literature 
for T > 100 K (Allison & Dalgarno 1969). In Fig. 9 we plot (5Tb as a function of time for recombinations starting at a = 30 
(solid lines), z = 20 (dashed lines), and z=15 (dotted). The initial temperature after reionization is taken to be Tk ~ K. 
The three curves in each case correspond to different values of the density contrast inside the bubble: 6=1 (higher curve 
), 6 = (middle curve), and 5 = —0.3 (lower curve). We have assumed uniform density inside the bubbles in solving the 
recombination equations. According to Fig. 9 a recombining bubble can yield a significant 21-cm emission of 5Th ~ 10 mK 
at (5 > for about 10*Yr for bubbles generated at 2; = 15 and 20. At « = 30 the emission lasts for a significant fraction of 
the Hubble time at that redshift. Electron-atom collisions are clearly important even for 5 = —0.3. Note that these bubbles 
can encompass large volumes of space so that the mean density in them is likely to be very close to 5 = 0. However, density 
variations inside a bubble will yield ST\, fluctuations at the level of differences between the curves shown in the figure. 



5 DISCUSSION AND CONCLUSIONS 

Recently, the cosmological 21-cm signal from the epoch of reionization has been the subject of intensive research activities. 
Twenty one cm maps are probably the only direct probe of the structure of high redshift reionization. In comparison, secondary 
CMB anisotropics will provide constraints on line of sight integrated quantities related to the ionized rather than neutral gas. 
Also, the statistical properties of 21-cm fluctuations contain valuable information on the physical nature of the ionizing sources 
and the amplitude of the underlying mass at high redshift. Further, it has become clear that the prospects for meeisuring 
the 21-cm signal are good in view of the various radio telescopes designed with this purpose in mind. Despite foreground 
contaminations (e.g. Oh & Mack 2003, Di Matteo, Ciardi & Miniati 2004, Morales & Hewitt 2004, Zaldarriaga, Furlanetto & 
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Figure 8. Top: The DBT fluctuations in a segment equivalent to L = 40 h^^Mpe in tlie frequency axis. All curves are for the same 
segment but as seen with different angular resolutions. Bottom: Correlation functions computed from the DBT as a function of frequency. 
In both panels, all curves for frequency resolution of Av = 100 KHz, and the solid, dotted, dashed, and dash-dotted lines correspond to 
angular resolutions of AO = 0, 0.4, 1.6, and 6.4 arcmin, respectively. 



Hernquist 2004), the collective data obtained by telescopes like the Low Frequency Array ^ (LOFAR), the Primeval Structure 
Telescope ^ (PAST), the Square Kilometer Array§, and the T-Rex instrument^ should, at the very least, help us discriminate 
among distinct reionization scenarios. 



t http://www.lofar.org 

•t http://astrophysics.phys.cmu.edu/ jbp 

§ http://www.skatelescope.org 

^ http://orion.physics.utoronto.ca/sasa/Download/Papers/poster/poster.jpg 
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Figure 9. The DBT as a function of time in recombining regions assumed to liave been fully ionized at redsliift ^ = 30 (solid curves), 
20 (dashed), and 15 (dotted). The ionizing source is switched off during recombination. The three curves at each redshift, respectively, 
correspond from top to bottom to mean density contrasts of 5 = 1, 0, and —0.3 inside the bubbles. 

In some reionization scenarios, collisions of free electrons with H I atoms can lead to a significant 21-cm cosmological 
signal from regions with density contrasts (1 + 5){x^j + y^/VHii^ — x^j)) > 20[(1 + z)/10]^^ and > 0.8, where the y^/y^^ 
is in the range 10 — 20 (Field 1958). We chose a simplified version of the X-ray pre-reionization scenario as an example for 
computing the fluctuations in the 21-cm brightness temperature. In this scenario decoupling by collisions is likely to play 
a major role, at least sufficiently far from the sources. We also demonstrated that collisions are also important during the 
recombination of gas in bubbles ionized by short lived sources. A recombining bubble should show in 21-cm emission (of 
JTb ~ 10 mK) for about lO^Yr. A thorough analysis of the collective emission from these bubbles is currently underway. 

In order to compute the fluctuations in the X-ray scenario, we have employed a simple semi-analytic model that allows 
us to assess the statistical properties of 21-cm brightness temperature fluctuations. The 21-cm fluctuations in this and other 
pre-reionization scenarios have been considered previously by several workers in the field (Bruscoli et al. 2000, Benson et al. 
2001, Liu ct al. 2001, Ihev et al. 2002, Ciardi & Madau 2003, Bharadwaj & Sk. Saiyad 2004, Ricotti, Ostriker & Gncdin 2004, 
Zaldarriaga, Furlanetto & Hernquist 2004) . In particular Ricotti, Ostriker & Griedin (2004) used numerical simulations of a 
1 h~^Mpc^ box to model various details of X-ray ionization scenarios. However, an assessment of the importance of collisions 
vs LyQ pumping as a function of position in the simtilation remains to be done. Further, the simulation box is still small and 
an analysis of the general statistical properties of the fluctuations must rely on semi-analytic models. 

We advocate here the the application of the Alcock-Paczyhski (AP) test (Alcock & Paczyriski 1979) on the correlations 
measured in three dimensional maps of 21-cm emission. These correlations are readily expressed in terms of angular separations 
and frequency intervals. In order to derive the correlations in terms of real space separations (measured for example in h~^Mpc) 
one needs to know the mean density parameter, Jim, and the cosmological constant Q.\. Constraints on these parameters can 
then be obtained by demanding isotropy of correlations in real space. A detailed analysis of the 21 cm AP test is presented 
elsewhere (Nusser 2004). 
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